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Abstract 

Patterned surfaces with large effective slip lengths, such as super-hydrophobic 
surfaces containing trapped gas bubbles, have the potential to reduce hydrody- 
namic drag. Based on lubrication theory, we analyze an approach of a hydrophilic 
disk to such a surface. The drag force is predicted analytically and formulated 
in terms of a correction function to the Reynolds equation, which is shown to be 
the harmonic mean of corrections expressed through effective slip lengths in the 
two principal (fastest and slowest) orthogonal directions. The reduction of drag is 
especially pronounced for a thin (compared to texture period) gap. It is not really 
sensitive to the pattern geometry, but depends strongly on the fraction of the gas 
phase and local slip length at the gas area. 

1 Introduction 

It is more than 100 years since Reynolds published his famous analysis of hydrody- 
namic lubrication due to a thin liquid film confined between two moving solids.— This 
theory is based on the simplification of the Navier-Stokes equations of continuum hy- 
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drodynamics by exploiting the special geometry of a film and no-slip boundary condi- 
tions. The utility of this theory still stands today and many extensions and applications 
of the original analysis maybe found in numerous publications. 

Classical solutions of creeping flow equations of the lubrication theory for a cir- 
cular disk of radius R moving with a velocity U towards a smooth wall, the so-called 
Reynolds problem, gives 1 

3 R^ oHU 
F R = ^nr]U—r, provided < 1. (1) 

when the gap H becomes small compared to R. Here r\ denotes the fluid viscosity. A 
consequence of this lubrication effect is that the close approach a disk to the wall, or 
its pulling away from it would take large time. This is the basis of the phenomenon 
of viscous adhesion, used in adhesives such as 'Scotch tape' or in the 'wringing' to- 
gether of smooth metal surfaces, 2 but may represent a very unfavorable scenario for 
other applications. An efficient strategy for reducing the drag force is to exploit hy- 
drodynamic slip, which can be generated at hydrophobic surfaces and is quantified by 
the slip length b (the distance within the solid at which the flow profile extrapolates 
to zero). 3 " 5 The near-field hydrodynamic interaction of a hydrophilic disk with such a 
hydrophobic surface (a situation which allows one to avoid a formation of a gas bridge 
in the gap 6 ) leads to a correction to the Reynolds force 7 

* F H + b 
f = — = — — (2) 
J F R H + 4b v ; 

Depending on the ratio b/H, the correction for slippage /* can turn to 1 (large com- 
pared to slip length distances) or 1/4 (small distances). Since for hydrophobic smooth 
and homogeneous surfaces b can be of the order of tens of nanometers ,— ~— but not 
much more, it is impossible to benefit of such a nanometric slip at separations 0(/im) 
and larger. 

Hydrophobicity can be significantly amplified by roughness, and extreme hydropho- 
bicity can be generated with well-controlled textures. 12 Such super-hydrophobic (SH) 
surfaces in the Cassie state, i.e. where the texture is filled with gas, can reduce friction 
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due to trapped nanobubbles, ' leading to a many-micron effective slip lengths.—— 
A mechanism for large slippage involves a lubricating gas layer of thickness e with 
viscosity r\ g much smaller than that of the liquid, the so-called 'gas cushion model' 18 



Taking into account that tj /rj g ~ 50, the variation of the SH texture height, <?, in the 
typical interval 0.1 — 10 jim gives b = 5 — 500 jim. The composite nature of the texture 
requires regions of very low slip (or no slip) in direct contact with the liquid, so the 
effective slip length of the surface, b e s, is smaller than b. Still, one can expect that a 
rational design of such a texture could dramatically reduce the hydrodynamic force at 
relatively large O(10/im) distances. 

Previous theoretical investigations of a flow past SH surfaces have addressed the 
questions of effective hydrodynamic^^ - — and electro-osmotic 2223 slippage. We are 
unaware of any previous work that has studied how the squeeze film drainage between 
surfaces would be modified by the occurrence of the effective slip. The main differ- 
ence from the simple model of a constant slip length, Eq. ©, used before for a smooth 
isotropic hydrophobic surfaces is that the effective slip length is itself not a character- 
istic of a heterogeneous (and, in general case, anisotropic) liquid/wall interface solely, 
but also depends on flow configuration, which in turn is determined by the smallest 
length scale of the problem, H, b, or roughness periodicity, L. Of course this imme- 
diately raises a difficulty: the decrease in H during the hydrodynamic interaction with 
SH surfaces would inevitably modify b e ff- 

In this paper, we explore what happens when a hydrophilic disc is driven towards a 
SH surface in the Cassie state. After describing the general theory in the following sec- 
tion, we present the results and discussion of the effect of SH slip in case of anisotropic 
and isotropic textures. Then follows a concluding section. 




(3) 



3 




Figure 1: (a) Sketch of a hydrophilic disk approaching a superhydrophobic surface, 
where a texture is represented in terms of patches of flow boundary conditions; (b) 
Examples of anisotropic (array of grooves) and isotropic (pillars) textures. 

2 General theory 
2.1 Model 

We consider a circular hydrophilic disk of radius R, which is parallel to and at a small 
distance H (<g; R) from a SH plane. Surfaces are immersed into viscous Newtonian 
liquid, and the pressure at the edge of the disk is atmospheric (p = po). The disk 
moves towards a plane with a constant speed U. This motion gives rise to an opposing 
force on the disk, which we aim to calculate. 

We examine an idealized SH surface in the Cassie state sketched in[T]where a liquid 
slab lies on top of the surface roughness. The liquid/gas interface is assumed to be flat 
with no meniscus curvature, so that the modeled SH surface appears as a perfectly 
smooth with a pattern of boundary conditions. The latter are taken as no-slip (b\ = 0) 
over solid/liquid areas and as partial slip (bi = b) over gas/liquid regions. We denote 
as 8 a the typical length scale of gas/liquid areas. The fraction of solid/liquid areas 
will be denoted 0i = [L— 8)/L, and of gas/liquid area <p2 = 1 — <pi — 8/L. Overall, the 
description of a SH surface we use is similar to those considered in Refs. 14 i 19 i 20 i 24 i 25 
In this idealization, some assumptions may have a possible influence on the friction 
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properties and, therefore, a hydrodynamic force. First, by assuming fiat interface, we 
have neglected an additional mechanism for a dissipation connected with the meniscus 
curvature. 26-28 Second, we ignore a possible transition towards impaled (Wenzel) state 
that can be provoked by additional pressure in the liquid phase. 29 ' 30 Both effects are 
expected to modify a hydrodynamic resistance force.—"— 

2.2 Governing equations 

We chose a Cartesian coordinate system with the origin at the SH surface, and z-axis 
directed towards the center of the hydrophilic disk. In our case of a thin liquid film 
and small Reynolds numbers (Re <C 1) they can be substantially simplified since the 
lateral component of the velocity field is large as compared with the normal component 
(v% ^> v,), and (V T v T ) <C dv T /dz- The Navier-Stokes equations are then reduced to 



where p is pressure, v T = v x e x + v y e y is the lateral velocity, and V T is the differential 
operator in plane (x,y), given by 

d d 

V T = — e x + — e v (5) 
ox ay 

The continuity equation takes the form: 



^r + (v T v T ) = o. (6) 

oz 



At z = H we have usual no-slip condition, while the boundary condition at z = 
reflects a tensorial hydrodynamic slip^ 



z = 0: (v T ) i = b iJ -^i,v z = 0; (7. 



z = H: ^ = 0,v z = -U. (8) 
Here and below we use the index form of representation for tensors and vectors: by 
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assuming summation over a repeated index, all indices (as like i, j, k, etc.) can be 
equal to either x or y. In particular, the first equation of © takes a form 

d 2 (v r )i 

n-^r 1 = Vip(*>y)> (9) 

where V; = (V T )„ 

2.3 Analysis 

Expression (|9} can be integrated twice over z, yielding the general solution for lateral 
velocity components (v T ),, i = x,y. In the classical Reynolds problem this gives two 
scalar constants, which should be determined via boundary conditions. Since in general 
case the SH texture is anisotropic, we have to assume the tensorial character of these 
integration constants, and find 

(vT)i=^(z 2 ^-AyZ-2?y). (10) 

Here 5y is the Kronecker delta (two-dimensional), Ay and By are constant tensors, that 
can be determined from conditions and ([8]). First, we find 

Then by substituting (fTTT) into ^} we get the tensorial relation 



B ik = b,jA jk . (12) 
Finally, by using condition © together with ( fT2l we get 

A ik H + b,jA jk =H 2 8 ik . (13) 

Equations (fT2l) and (13[ determine unknown constants in the expression for a tangential 
velocity, Eq. (TTOb . 
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To simplify further analysis, we now align basis vectors with principal directions 
of the slip length tensor {bjj} 



b e ff : 











bis 



(14) 



where the eigenvalues £r ff an d b^ s are effective slip lengths in the fastest and slowest 
directions, correspondingly. These values can be related to the components of the 
effective channel permeability tensor, that determines the average fluid flux across the 
channel's cross-section, and, thus, depend on the gap thickness H. 

Now we can explicitly calculate components of {A,y} and {Bjj} and conclude that 
their principal directions coincide with those of {bij} tensor 



H 1 



H + b eff 



"j A>£y Aw 0; Ayy 



H 1 



(15) 



Byoc — b ef{ 



H 1 



H + b 



; Bxy — B yx — 0; B yy = b eS 



H 2 



eff 



By integrating the continuity equation ^ 



(16) 



li 

U = J(V T v z )dz, 



(17) 



we obtain the expression for the relative speed of surfaces 



(18) 



which represents a partial differential equation for pressure: 



^ d 2 p „ d 2 p 

r — 4- C — 

Lx dx 2 + y dy 2 



-u, 



(19) 



where 



H 3 H + 4b 
12TJ H + b 



cff 



cff 



(20) 
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and 



_ H 3 H + 4b^ 
* 127] H + b^ 



(21) 



The exact solution of this partial differential equation satisfying the boundary con- 
dition for pressure at the edge of disk is 



P=P0 + 



U (R 2 



2 (C X + Cy) 



2 i 2 



(22) 



The hydrodynamic resistance force F acting on the hydrophilic disk of radius R is 
opposite to the force exerted by the SH surface. We remark and stress that although the 
anisotropy of a texture leads to the reduction of the physical symmetry of the whole 
system, the resulting force is still directed along the axis of highest symmetry, e z . Thus, 
lateral force components F x = F y vanish due to the presence of mirror planes parallel to 
the z-axis. The only remaining component of the drag force is denoted below as F, and 
in the first-order approximation may be evaluated as the integral over the disk's surface 



F =2% J \ p-pa-2T]-^- )rdr 

o 



(23) 



However, it first-order approximation we may omit the last term in the integrand, and 
obtain— 

„ 3 nriUR 4 



2 H 3 

where the correction for an effective slip is 



"/eff; 



(24) 



ft 



eff 



Fr 



H + Ab\{H) H + 4b^{H 



H + bl s (H) H + b^(H) 



-l 



(25) 



Thus the effective correction for a SH slip is the harmonic mean of corrections ex- 
pressed through effective slip lengths in two principal directions, 



1 

/eff 



Hi 



i 



eff 



/eff , 



(26) 



In case of an isotropic textures, all directions are equivalent with b\ s — b^ s = b e ff, so 
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Figure 2: Special textures arising in the theory: (a) stripes, which attain the Wiener 
bounds of maximal and minimal effective slip, if oriented parallel or perpendicular to 
the pressure gradient, respectively; (b) the Hashin-Shtrikman fractal pattern of nested 
circles, which attains the maximal/minimal slip among all isotropic textures (patched 
should fill up the whole space, but their number is limited here for clarity); and (c) the 
Schulgasser and (d) chessboard textures, whose effective slip follows from the phase- 
interchange theorem. 



we get 



/eff — 



F _ H + b eS (H) 



(27) 



/efr F R H + 4b eS {H) 

Note a similarity to Eq. (0. The only difference is that slip length is effective and 
dependent on H. Obviously, the case b\ n = b^ s = corresponds to f* s = 1 and gives 
Eq. O 



3 Results and discussion 

In the preceding section, we derived a general expression for f* s , which relates it to the 
effective slip length of the SH wall and the gap. In order to quantify the reduction of a 
drag force due to a presence of SH wall, this expression will now be examined for some 
specific anisotropic and isotropic textures, where analytical or numerical expressions 
for bh^ have been obtained. 



3.1 Anisotropic textures 

To highlight effects of anisotropy and to explore the effect of gap on the force, we now 
focus on flat, periodic, striped super-hydrophobic surface (sketched in|2a)), where the 
local (scalar) slip length b varies only in one direction. Such surfaces have been already 
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used for reduction in pressure-driven flows 15 and enhancement of mixing. 36 

The problem of flow past stripes has also been examined theoretically. Effective 
slip lengths (H) in case of an arbitrary channel thickness can be calculated semi- 
analytically following approach, 28 and the details will be published elsewhere. In the 
case of thick channel (H L) the effective hydrodynamic slip reads 20 



y eff P 



L 

K 



In 


( x<h \ 






L 

1 + — In 
nb 


sec 
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+ tan 


( 





(28) 



In 

L 


sec 
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H In 
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sec 


( K<h\ 


+ tan 


( 





(29) 



Flow in a large channel does not depend on H, and is controlled by the ratio of the local 
slip length b to texture period L. At b/L ^> 1 expressions (|28]|-(|291 turn to 



y eff 



-In 



% <p 2 



and b£ e - 



y cff 



(30) 



suggested earlier for a perfect (b = «>) local slip. 21 ' 23 i 25 ' 28 ' 37 At b/L < 1 Eqs. (l28T>-d29b 
predict a simple surface average isotropic effective slip 



J eff 



b<p 2 . 



(31) 



In the case of thin channels (H <C £) striped surfaces were shown to provide rigor- 
ous upper and lower Wiener bounds on the effective slip over all possible two-phase 



pattern: 



19 



|| bHfa 



bH<p2 



n H + 4b(pi 

At b/H ^> 1 these give truly tensorial anisotropic effective slip 



(32) 



y eff 



H-, b sB - 



"eff 

4 ' 



(33) 
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Figure 3: Eigenvalues bn/H (solid curve) and b±/H (dashed curve) of the normalized 
slip length tensor for stick-slip stripes of period L with local slip length at liquid-gas 
interface b/L = 1 and slipping area fraction (j>2 = 0.5 as a function of H/L. 

but at b/H <C 1 it leads to Eq. (|3T1 i. The above formulae well illustrate the fact that ef- 
fective boundary conditions are controlled by the smallest length scale of the problem. 

Thus, the effective slip, Z? e ff, is large at H/L 1, and decreases when H/L is 
small. However, according to Eq. (f25l to reduce a drag force we need to maximize the 
ratio b c ff/H, but not the absolute values of effective slip itself. The computed results 
for bh^/H presented in |3]for b/L — 1 and (j>2 = 0.5 show that these values become 
discernible when H/L = 0(1) and smaller, by giving asymptotic values predicted by 
Eqs. ( 1331 ). They however vanish at large H/L, which is the consequence of the fact 
that according Eqs. (l28ll-(|29i> the effective slip length, in this geometry, is mainly fixed 
by the texture period, L, so that by slightly modifying the results of~ we get b^g/H °< 
— (L/H)ln(l — (fe). This suggests that for feasible (02 below 0.99) surfaces a significant 
reduction of hydrodynamic drag would be possible to obtain only at a thin gap limit. 

This is illustrated in [4] where other calculations of b e ff/H, made for several 02 an d 
b/L, were used to compute the correction for effective slip, /* ff , as a function of the 
gap. Indeed, at large distances all curves converge to /* ff =1. In other words, the 
drag force is the same that it would be in case of a hydrophilic surface, F — Fr. This 
conclusion can be derived directly from Eq. d25b and is valid for any, however large, 
b. Results presented in [4] show that for a thin gap the correction for effective slip has a 
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Figure 4: The correction for effective slip, /* ff , versus dimensionless gap width H/L for 
striped superhydrophobic surface. Solid curves correspond to local slip length b/L = 
10 (from top to bottom 02 = 0.2, 0.5 and 0.9), dashed curves - to b/L = 0.1 (from top 
to bottom 2 = 0.2 and 0.5), dash-dotted curve - to b/L = 0.01 and 2 = 0.5. 

tendency to decrease with 02. Substitution of Eqs. (f32t into Eq. d2"5b allows to quantify 
this important result 



2(H + 4b-3b<j> 2 )(H + b) 

/eff 2H 2 + lQbH + Eb 2 + 9b 2ty 2 _ 9b 2fy2 



(34) 



In case of a small local slip, b/H <C 1, we derive 



eff 



1-3 — 02 



(35) 



However, a more important limit that would represent a minimal possible, but feasible, 
value of f* s for a striped texture of a given 02 can be attained in case of a large local 
slip, b/H^>l 



eff - 



2(4-302) 

8 + 90 2 - 90 2 2 



(36) 



This expression, in particular, shows that f* s varies in the interval from 1 to 1/4 for 02 
between and 1, which is in agreement with initial expectations. 
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Figure 5: Effective slip length, b e ^/H, versus b/H [for <p2 = 0.5] in a thin gap limit, 
H <^L. Superhydrophobic surfaces are: anisotropic stick-slip stripes attaining Wiener 
bounds (dashed curves), isotropic textures attaining Hashin-Shtrickman bounds (solid 
curves) and satisfying the phase-interchange theorem (dash-dotted curve). 

3.2 Isotropic textures 

Consider now isotropic structures, without a preferred direction. Textures such as 
arrays of pillars (posts) or hollows represent a very important experimental geome- 
try. V 3 i 38 No exact analytical or semi-analytical solution of the Stokes equations has 
been performed up to now for isotropic patterns for a gap of arbitrary thickness or even 
in the limit of thick channel. Based on numerical results^ 4 - we conclude that all these 
textures provide effective slip confined between expected for transverse and longitudi- 
nal stripes except as in the limit of vanishing solid area (02 ^ 1)- In the latter case, 
scaling arguments 5 suggested for patterns of individual pillars, where the largest effec- 
tive slip is expected, b e s/H °c (L/H)/ (n^/l — fa). Obviously, with the realistic <p2 this 
cannot change the above conclusion made for anisotropic surfaces, that a significant re- 
duction of hydrodynamic drag would be possible only for a thin gap. Therefore, below 
we focus on a thin gap situation, by trying to highlight the effect of texture geometry. 

Rigorous (Hashin-Shtrikman) upper and lower bounds on an effective slip length 
for arbitrary isotropic textures, given only the area fraction and local slip lengths, where 
calculated analytically in. 19 i 39 These bounds can be attained by the special Hashin- 
Shtrikman fractal pattern 40 sketched in|2jb). For one bound, space is filled by disks of 
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all sizes, each containing a circular core of one component and a thick ring of the other 
(with proportions set by the concentration), and switching the components gives the 
other bound. Fractal geometry is not necessary, however, since periodic honeycomb- 
like structures can also attain the bounds.— For a situation considered here the corre- 
sponding upper bound for the effective slip length in case of isotropic surfaces can be 
deduced from the general resul t 19 > 39 as 



eff H(2H + 5b)+bfa(5H + 8b) K ' 



The lower bound is then 



2H + 5o0i 

At large b/H these give for upper and lower bounds 

5i/0 2 2H<j>2 
b ^=~^T ' and ^ff=-z— , (39) 
80i 50i 

correspondingly, i.e. similarly to anisotropic stripes (cf. Eq. [33l >. b e s/H scales as 
°< 02/ 01 ■ The bounds for b e s /H are plotted versus H/b in [5] Also included are results 
for Wiener bounds, Eqs. (l32l l. plotted in a same way. Finally, for completeness we add 
phase interchange results, which in the particular case of a medium which is invariant 
by a 7C/2 rotation followed by a phase interchange, gives 

i.—j—.-H, (40, 

Obviously, 0i = 02 = 0.5 for such a medium. Classical examples of such an isotropic 
texture are the Schulgasser proposal, or a family of chessboards, examples are shown 
in Etc) and (d). If b/H 3> 1 for textures that follow phase interchange theorem we 
simply get 

£eff=f (41) 

[5] shows that Hashin-Shtrikman bounds are relatively close and confined between 
Wiener ones. It can also be seen that all curves behave similarly, by vanishing at large 
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Figure 6: The correction for effective slip, /* ff , versus b/H [for 02 = 0.5] in a thin gap 
limit, H <C L. Superhydrophobic surfaces are: anisotropic stick-slip stripes attaining 
Wiener bounds (dashed curve), isotropic textures attaining Hashin-Shtrickman bounds 
(solid curves) and satisfying the phase-interchange theorem (dash-dotted). 



H/b. This is a consequence of the fact that in this limit the effective slip coincides with 
the average, fo e ff — fob (cf Eq. (|3T1>). At small H/b all plotted curves give a plateau de- 
scribed by Eqs. ( 1331 ). d39l l. and (HTt . depending on the texture. Its height is controlled 
solely by 02/01 an d texture type. 

Using Eq. d27l ) together with Eqs. ( f37l >. ( f38l > we obtain the corresponding lower 



eff 



(H + b)(8b-3fob + 2H) 
{H + 4b)(2b + 3<j) 2 b + 2H) 



and upper 



2H + 5b-3<j) 2 b 



eff - 



(42) 



(43) 



2H + 5b + 3fob 

bounds for a correction for effective slip. We remark here, that the lower bound for 
/* ff corresponds to the upper bound for b e s, and vice versa. Similarly, by combining 
Eq. d27b together with phase interchange results, Eqs. (1401 1. we derive 



eff 



H + b 
H + 4-b 



(44) 



The Hashin-Strickman bounds for f* s are plotted in |6] The results for textures 
satisfying phase interchange theorems are also included in [6] and confined between 



15 



Hashin-Strickman bounds as predicted by the theory. To examine the significance of 
isotropy/anisotropy more closely, the short-distance region of one of the curve from[4] 
is reproduced in [6] in the corresponding coordinates. It turns out that the results for 
stripes are also confined between Hashin-Strickman bounds for /* ff , so that in general 
case anisotropy would not help to reduce/enhance a drag force. We stress however, that 
the Hashin-Strickman bounds are fairly close, so the theory provides a good sense of 
the possible f* s of any isotropic or even anisotropic texture. The data presented in [6] 
show that at very large distances, the resistance to approach flow is the same as it would 
be in the Reynolds problem, with no slippage on both surfaces. A straightforward 
calculation shows that at small b/H the useful approximation for f* s for all textures 
would be Eq. ( f35b - This universal behavior is confirmed by coincidence of all curves 
presented in [6] in this limit. If the gap is much smaller than local slip length at the 
gas area, the correction for effective slip becomes smaller and turns asymptotically to 
constant values. For the Hashin-Strickman bounds these can be evaluated as 



8-3<fe 5 - 30 2 

4(2 + 302)' eff " 5 + 3<fc 



/rff- I ~ , /eff ^^Tf (45) 



Correspondingly, for a Schulgasser (or chessboard) textures /* ff ~ 1/2 in this limit, 
which is well seen in|6] 

These results suggest that the key parameter determining reduction of drag is the 
area fraction of gas, 02, in contact with the liquid. This is illustrated in|7l where (using a 
relatively large b/H) Hashin-Shtricknann bounds for f* s are plotted versus the liquid- 
gas area fraction 02- If this is very small (or 0i — > 1) for all textures, the correction 
for slip tends to its absolute maximum, f* s = 1. In the most interesting limit, 02 —> 1, 
we can achieve the minimum possible value of correction for effective slip, f* s —1/4, 
provided b/H is large enough. We also notice that at small 02 the results for stripes 
coincides with the lower Hashin-Shtricknann bound for /* ff . In contrast, while for large 
02 they reduce a drag as it would be in the upper Hashin-Shtricknann bound for /* ff . 
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Figure 7: The correction for effective slip, f* s , versus 02 [f° r b/H = 15] in a thin gap 
limit, H <C L. Superhydrophobic surfaces are: anisotropic stick-slip stripes attaining 
Wiener bounds (dashed curve), isotropic textures attaining Hashin-Shtrickman bounds 
(solid curve for the upper bound and dash-dotted for the lower one ). The value for the 
chessboard or the isotropic Schulgasser structure is also shown as a circle. 

4 Conclusion 

We have analyzed the squeeze-film drainage of a liquid confined between a hydrophilic 
disk and patterned SH surface of non-uniform local slip length, and have obtained gen- 
eral solutions to arbitrary gap, and slip variation. We have shown that the decrease 
in the hydrodynamic force in the presence of a patterned slipping surface can be de- 
scribed in terms of a correction for slippage to the Reynolds formula, formulated as a 
function of the slip lengths in the fastest and slowest direction. Provided the separation 
is small compared to texture period, this correction to slippage becomes small as com- 
pared with unity, providing the significant decrease in hydrodynamic drag. We have 
concluded that in all situations, to achieve a large reduction of a drag force optimizing 
the pattern geometry is not nearly as important as to maximizing local slip at the gas 
area and the fraction of the gas phase. 
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